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Abstract 

Small objects can swim by generating around them 
fields or gradients which in turn induce fluid mo- 
tion past their surface by phoretic surface effects. 
We quantify for arbitrary swimmer shapes and sur- 
face patterns, how efficient swimming requires both 
surface "activity" to generate the fields, and surface 
"phoretic mobility." We show in particular that (i) 
swimming requires symmetry breaking in either or 
both of the patterns of "activity" and "mobility," and 
(ii) for a given geometrical shape and surface pattern, 
the swimming velocity is size- independent. In addi- 
tion, for given available surface properties, our calcu- 
lation framework provides a guide for optimizing the 
design of swimmers. 

Introduction. In the current miniaturization race 
towards small motors and engines, a rapidly expand- 
ing subdomain is the quest for autonomous swim- 
mers, able to move in fluids which appear very vis- 
cous given the small length scales (low Reynolds num- 
ber). Robotic microswimmcrs that generate surface 
distortions is an avenue (e.g. by mimicking sperms 
[1]), but it seems equally interesting to try to take 
advantage of physical phenomena that become pre- 
dominant at small scales. Intcrfacial "phoretic" ef- 
fects (electrophoresis, thermophoresis, diffusiophorc- 
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sis, [2]) are from this standpoint a natural avenue 
given the increased surface to volume ratio. 

Very recently, in line with earlier suggestions of 
self-electrophoresis of objects that would be able to 
generate an electric field around them [3J [U d] and 
our recent proposal of self-diffusiophoresis [5], many 
experimental reports have appeared of heterogeneous 
objects (e.g. rods of hundreds of nanometers) swim- 
ming using different catalytic or reactive chemistry 
at their different ends [HI [HI El IE] • Although many 
mechanisms may contribute to the observed motion 
[6], phoretic mechanisms have been shown to be es- 
sential for some systems [TO] . 

In this Letter we consider generally self-phoretic 
motion, where the swimmer, generates through its 
surface activity gradients of (at least) a quantity 
(concentration of a dissolved species, electric po- 
tential, temperature), which in turn induce motion 
through classical interfacial phoretic processes. Be- 
cause it is now possible to fabricate such objects with 
controlled shape and patterns of surface properties 
El [11] , we discuss how these design parameters af- 
fect the swimmer's velocity. A general framework is 
proposed to compute this velocity for arbitrary swim- 
mer size, shape and surface pattern. We then work 
out explicitly a few examples for spheres and rods, 
geometries amenable to analytical calculations, and 
corresponding to existing means of fabrication [51 [7] . 
From these examples we extract generic rules as to 
the role of the pattern symmetry, of the swimmer size 
and shape, before concluding with a few remarks. 
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General formalism. We start with a reminder of 
what phoretic mechanisms are [2]. Consider e.g. a 
colloid in a solution where the concentration of a so- 
lute c exhibits a gradient. The colloid is then set into 
motion due to the induced imbalance of osmotic ef- 
fects within the solid/liquid intcrfacial structure at 
the surface of the object. When the thickness of the 
interfacial layer is thin compared to the object, the 
resulting flow is most conveniently described by an 
effective slip velocity of the liquid past the solid, pro- 
portional to the local gradient of c [5] : 

v s (r s )=/i(r s )(I-nn)'Vc(r s ) (1) 

where n is the local normal to the surface, I is the 
identity tensor. /ti(r s ) is the local "surface phoretic 
mobility" , positive or negative depending on specifics 
of the solute/surface interactions [2]. 

This distribution of slip on the object's surface 
induces a net drift velocity V, which can be com- 
puted using the reciprocal theorem for low Reynolds- 
number hydrodynamics [12j . Its components in a ref- 
erential (ej) are: 

V • e, t = - JJ dr s n • cr; • v s (2) 

where er^ is the hydrodynamic stress tensor at the 
surface of an object of similar shape dragged by an 
applied unit force exactly equal to ej, in the absence 
of slip. 

Similar principles describe the motion of the ob- 
ject in a gradient of electric potential (electrophore- 
sis), or in a temperature gradient (thermophoresis). 
The above formalism can be applied upon simple re- 
placement of c by the potential <f> e i or by the tem- 
perature T. For electrophoresis /i is up to a con- 
stant the so-called zeta-potential [2|. For all these 
situations, achieving large values of fi requires en- 
gineering of the surface in relation to the solvent 
(e.g. nanometer scale roughness and solvophibicity 
do matter pi [14]). 

We now focus on objects that generate the gradi- 
ents themselves, e.g. using the chemical free-energy 
available in the solvent [6l [7]. Returning to self- 
diffusiophorcsis, the simplest situation corresponds 
to a solute generated by active sources and sinks on 



the surface of the object, and diffusing in the bulk. 
At steady-state, in the reference frame of the object, 
and neglecting distortions induced by the flow (small 
Peclet number), the solute concentration in the liquid 
is given by: 

L>V 2 c = 0, (3) 

-Dn- Vc(r s ) = a(r s ), (4) 

where D is the diffusion coefficient, and (x(r s ) mea- 
sures the "surface activity" at position r s on the sur- 
face, i.e. the generation or consumption of solute by 
a chemical reaction. In general, describing this pro- 
cess involves additional coupled transport problems 
for other species involved in the surface reactions. 
For simplicity, we nevertheless proceed with a fixed 
local value of a. This should well describe situations 
where the solvent contains an excess of the reactants 
necessary to produce and destroy the solute of inter- 
est, so that the density of active sites on the surface 
is what limits the fluxes, a is then this density times 
the rate of the chemical reaction per site [5] . 

For self-electrophoresis, similar formulae hold upon 
replacing c by </> e ; and D by the solution electrical 
conductivity, a is then the electrical current injected 
by the surface. Existing realizations [9l [10], corre- 
spond to no net current generation JJ S dr s a(r s ) = 0, 
with the ions produced on one side of the swimmer 
and consumed on the other one, while electrons are 
transported through the body of the swimmer. Sim- 
ilarly, for thermophoresis, beyond replacing c by T 
and D by the thermal conductivity, the boundary 
condition (4) must be adjusted to account for heat 
transport through the object. 

The general formal description (1-4) predicts a 
swimming velocity V ~ a/x/D, linear in the surface 
properties a and n, and inversely proportional to the 
medium conductivity D. More interestingly, once the 
pattern of surface properties and the object shape arc 
given, V is independent of the size R of the swim- 
mer. This central point emphasizes the robustness of 
this swimming strategy against downsizing. In con- 
trast, an external body force on the object scaling as 
F ~ R 3 (dielectrophoresis, magnetophoresis) leads in 
this viscous realm to velocities ~ R 2 , which decrease 
rapidly with length scale. Further, our approach (1- 
4) allows us to relate quantitatively the velocity V to 
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the design of the swimmer (shape, surface patterns 
of a and fi), as we show analytically for spheres and 
rods. 

Spheres. For a sphere of radius R, as ner; = 
(l/47ri? 2 )e;, equation (2) becomes 



V 



47 ~R2 jj dVs v(r s )(L - nn) • Vc(r s ), 



(5) 



We focus further on azimuthally symmetric patterns, 
with surface quantities depending only on the "lat- 
itude" angle 9 with the polar axis e 2 . We ex- 
pand in Legcndrc polynomials. The surface activity 
a (Q) = Sfco &iPe(cos 9), generates a concentration 
variation given by Eqs. © and ([J|: 



c(r,9) 



R 



E 



d i + i 

1=0 



R 



Pi(cos9), (6) 



where Coo is the concentration at infinity. For a sur- 
face mobility fi(6) = /i^P^(cos#), the corresponding- 
gradient yields a swimming velocity (5): 
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(7) 

As a first check, we recover [5] that a single ac- 
tive site at the pole a = t~ 1 62d(*s) of a sphere of 
uniform mobility /j, yields V — ilT ^i Dr ■ By linear- 
ity, an emitter at the back pole with rate l/r e and a 
consumer at the front pole with rate l/r c , results in 

V = tL ( -L 4- -L\ 

v TttWd \t c ' t c j- 

To access useful general principles we now focus on 
three simple designs (Figure 1), within reach of cur- 
rent fabrication techniques [llj . The first is the Janus 
sphere, with two hemispheres covered with a different 
enzymatic or catalytic material. Generically, their 
mobilities /i also differ, so we take: 



{a{6)^{6)) = 



and find analytically 

v = i V 

8 D K 




< 9 < 5, 



< 8 < 7T, 



(8) 



(9) 




Figure 1: Three spherical swimmers described in the 
text. The flux of c (opposite to the gradient) is de- 
picted with the dashed arrows, while the resulting slip 
velocity of the liquid along the sphere is described by 
open arrows. The net swimming direction is given by 
the thick dark arrow, (a) a Janus swimmer, with 
a homogeneous /i and a broken top-down symme- 
try in a: one hemisphere produces c while the other 
consumes it. (b) a Saturn swimmer with an equa- 
torial belt producing c, and a broken symmetry in 
the phoretic mobility between the two hemispheres 
(here even opposite signs), which induce swimming 
in the same direction from the (opposite) equator 
to poles gradients, (c) a three-slice design in with 
polar slices efficient in producing and consuming c, 
while the equatorial slice is chosen for its large mo- 
bility. Properly choosing 8q permits maximization of 
the swimming velocity (Fig. 2) 

This shows that symmetry breaking in the chemical 
activity (a_ ^ a+) is essential to propulsion, and 
that the swimming velocity is larger if the two mo- 
bilities arc of the same sign (which dictates the swim- 
ming direction). 

The second design is a "Saturn" particle, where 
the surface activity is concentrated around the equa- 
tor, and where symmetry is now broken by choos- 
ing hemispheres of different mobilities. For computa- 
tional convenience we solve the very similar problem 
defined by: 



a(9) = o>q(1 — cos 2 9), fi{6) = /j,q cos 9, 



and find 



V 



4 n a 



(10) 



(11) 



45 D 

So symmetry breaking in mobility alone also leads 
to swimming: the activity at the equator generates 
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gradients of c towards the poles that generate slip 
velocities in the same direction along z as the sign of 
/i is opposite on the two hemispheres (Figure lb). 

The third design is that of a three-slice sphere, 
that illustrates typical considerations that arise when 
trying to increase by design the swimming velocity. 
From the scaling analysis above one should obviously 
aim for surfaces with both large values of a and /i. 
However, in many cases these goals will not be met 
by a single surface, so one may instead look for ways 
to combine in an efficient design surfaces with large a 
and surfaces with large [i. As an example, let us try 
to upgrade the Janus design of Figure la by adding a 
passive belt of large phoretic mobility /io, which leads 
to the swimmer of Figure lc: 



\V\/[a^/D) 



(<*+,//+), 0<6< 



+ ■ 



(0,Mo), 
(a-)M-)) 



< 6 < 7T - 



< 6 < 7T, 



(12) 

To reach the largest swimming velocity for given 
surface properties, one can adjust the values of 6 + 
and Q- in a compromise between increasing the mo- 
bility by enlarging the belt without shrinking too 
much the poles that generate the gradient. Figure 
2 shows results for the model case a + = —at- = a, 
/i + = /!_ = fi of same sign than /^o, 9 + = 6- = 9q\ 
for /io > ^ there is a value of 9q that maximizes the 
speed, which is always larger than 7r/3, i.e. no matter 
how large /i , the " optimal" belt should not be made 
larger than a fixed value. The relative flatness of the 
bottom curve in Figure 2 shows a moderate depen- 
dence on #o around this value, so that the swimming 
velocity should be rather robust to small imprecisions 
in the fabrication. 



Thin rods. Because of the many experimental re- 
alizations of such systems [H [7] , we discuss the case 
of long cylindrical rods ( length 2L and radius b <C L) 
with surface properties varying only along their axis 
—L < z < L, i.e. fi = n(z),a = a(z). In the limit 
L > d slender body theory can be used, and remark- 
ably equation |(2J) yields a simple form for the swim- 




Figure 2: Reduced swimming velocity V/(a^o/D) of 
the three-slice swimmer of Figure lc, as a function of 
"polar" angle 8o for various values of (J,/fJ,o: 0, 1, 4 
(from bottom to top). The optimal value at low \ij \iq 
reflects the compromise between large poles to have 
strong chemical gradients, and a large equatorial belt 
for these gradients to translate into large flow motion. 
Inset: optimal angle 9q against n/jJ,o, ■ 

ming velocity along z: 

V ^ J ^dz fi(z)d z c s , (13) 

In the same limit, the surface concentration c s is re- 
lated to the surface activity a(z) by [6] 

b f L , a(z') 

Cs{z) ^2Dj_ L dz V^-y (14) 

These two equations give the velocity of a rod for any 
"bar-coded" surface properties a and fi. For a Janus 
rod 

{(a_,/x_), — L < z < 0, 
(15) 
< z < L, 

we find 

which illustrates a few important points: (i) again 
swimming is obtained by symmetry breaking in the 
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chemical activity or/and in the phoretic mobility, (ii) 
swimming with uniform activity (a ^constant) is 
possible - in contrast to the special case of spheres 
for which this leads to no tangential gradients of c 
-, (iii) the velocity scale is decreased by a factor of 
~ b/L when compared to a sphere of radius ~ L, 
because the surface generated flux of c is reduced 
by a similar factor, and (iv) V is up to logarithms 
proportional to 1/L at fixed rod diameter, but scale 
independent if length and diameter are varied in the 
same proportion. 

Discussion. The swimming velocity V of a 
phoretic swimmer, is independent of its size and 
scales as V ~ afi/ D, where a and /x are some aver- 
ages quantifying the surface activity generating the 
gradient and the surface phoretic mobility. The exact 
value of V depends on the shape and surface pattern 
in a computable way. For example, spheres swim 
faster than rods of same longest dimension, and im- 
provement can be obtained by properly combining 
various types of surfaces. 

- This type of formalism holds for self- 
diffusiophoresis (with a the surface reaction density, 
\x the diffusiophorctic mobility, and D the solute 
diffusion coefficient), for self-electrophoresis (with a 
the electrical current density injected in the solution, 
\x the electrophorctic mobility, and D the solution 
electrical conductivity), and for self-thcrmophoresis 
(with a the surface heat flux injected in the solution, 
\i the thcrmophorctic mobility, and D the solution 
thermal conductivity). As for motion driven by ex- 
ternal gradients, the largest velocities are expected 
for highly charged surfaces at low ionic strengths, 
which yield large mobilities /i because of the thick 
interfacial electrical double-layer. 

- The present line of reasoning can be adapted to 
more accurate descriptions of surface chemistry for 
specific experimental conditions. What is required in 
a given situation is an appropriate scheme for com- 
puting a(r s ). For example if c is produced at the 
surface from some chemical at concentration cj^ el far 
away from the object, various regimes show up. We 
have focused on kinetics limited by the number of 
catalytic sites on the surface, but for low values of 



c fuei t ne reaction becomes diffusion limited so that 
a depends on the global geometry. 

- Due to thermal rotational diffusion, if no exter- 
nal field orients the objects, the direction of mo- 
tion persists only for a time r ~ l/D rot , where 
D rot ~ ksT/riB? and R is the largest dimension 
of the object. Note however that the field c typ- 
ically adjusts much faster to a change of orienta- 
tion (R 2 /D <C l/D rot ), so our steady-state esti- 
mates still hold for the instantaneous velocity. A 
simple Langevin dynamics with a diffusion coefficient 
Dtrans ~ k B T j r\R and drift at V along a stochasti- 
cally changing direction (D rot = 1/r) leads to diffu- 
sive behaviour at long times (f>r) with an activity- 
enhanced diffusion coefficient D e ff = Dt ra ns + ^V 2 T. 
This enhancement may be significant and measurable 
for not too small swimmers, i.e. R > {k B T /rfV) 1 / 2 , 
of order ~ a few 100 nm for V ~ 10 ^m.s -1 in water. 

For nano-rods of radius b fixed due to fabrication 
constraints [6] , the scaling of the previous quantities 
become (using the shorthand Vb = a/fiD and forget- 
ting logarithms) V ~ V (b/L), r ~ ijL 3 /k B T, and 
an activity-induced enhancement of the bare trans- 
lational coefficient Dtrans ~ kBT/rjL by a quantity 
~ (r)V^/k B T)b 2 L, significant if L/b > (k B T/r)V b 2 ). 

- Beyond the "linear" swimmers described here, 
our approach can also describe active "spinners," 
which rotate using the same mechanisms [6l [7] . Their 
angular velocity can be computed by replacing (4) by 
the formula for angular swimming speed given in Rcf. 
|12j . and scales as f2 ~ a^i/DR, decreasing with ob- 
ject size for a given shape and geometrical surface 
pattern. 

- There is a growing interest in the collective be- 
haviour of swimming organisms (see e.g. [15j and ref- 
erences therein). Phoretic swimmers are special with 
regards to these questions as they interact through 
both their hydrodynamic fields (decaying as the in- 
verse of the cube of their relative distance r -3 ) [2] 
and the c-ficld gradients they generate (which decay 
as r~ 2 for a net production/absorption of c per ob- 
ject, and as r -3 otherwise). 

In conclusion, we have provided generic considera- 
tions for the design of small phoretic swimmers (size 
(in)dependence of the swimming velocity, necessary 
symmetry breaking in the surface pattern), as well 
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as quantitative procedures to estimate the influence 
of their velocity. Beyond these considerations, fast 
motion of course relies on clever surface engineer- 
ing to obtain large surface "activity" and "phoretic 
mobility." We hope that our calculations will stim- 
ulate experimental studies of swimmers with given 
surface chemistries but various surface patterns, e.g. 
nanorods U\ and spheres [TTJ [TB] , keeping in mind 
that other effects (e.g. bubble generation [7]) may 
compete with the ones described here. 
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